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A general tool for description of open quantum systems is given by the formalism of quantum 
operations. Most important of them are trace-preserving maps also known as quantum channels. 
We discuss those conditions on quantum channels that the Jarzynski equality and related fluctuation 
theorems hold. It is essential that the representing quantum channel be unital. Under the mentioned 
condition, we first derive the corresponding Jarzynski equality. For bistochastic map and its adjoint, 
we further formulate a theorem of Tasaki-Crooks type. In the context of unital channels, some notes 
on heat transfer between two quantum systems are given. We also consider a finite quantum system 
operated by an external agent with a feedback control. When unital channels are applied at the 
'_^»T ■ first stage and, for a mutual- information form, at the further ones, we obtain quantum Jarzynski- 

Sagawa-Ueda relations. These are extension of the previously given results to unital quantum 
^— ^ , operations. 
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0\ ; I. INTRODUCTION 

Evolution of complex systems away from thermal equilibrium is one of fundamental issues of statistical physics 
[ltij . The growing interest to mechanics and thermodynamic of systems at the nanoscale has also stimulated a more 
detailed study of statistical fluctuations 15.]. In addition, the role of computer simulations in analyzing nonequilibrium 
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processes will certainly increase. For instance, many efficient schemes for simulation of complicated biomolecules 
and colloidal particles were developed 35 1. Due to Jarzynski [ll], EH; important advances have been achieved for 
q. thermodynamic systems driven out from the equilibrium by external forces [TBI l37j . If these forces are varied in line 
—t ' with a specified protocol, then some exact relations can be derived. The first of exact nonequilibrium relations is now 
referred to as the Jarzynski equality. Results of such a kind are significant in own rights as well as for extending the 
scope of computer simulations (26|. New relations have been tested in experiments (l9l lifi |42| . Various aspects of 
Jarzynski's equality and related fluctuation theorems have been addressed in the papers [22|, QM, [33, EH . The issue is 
still the subject of active research 0, [T3, [H, [33[ . Further details and bibliography can be found in the review 15 1. 
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We consider a thermally insulated system, which is acted upon by a time-dependent external field. For any quasi- 
static process, the work performed on the system concurs with difference between the final and initial free energies. 
For a nonequilibrium process, averaged total work will exceed this difference. The claim directly follows from the 
second law of thermodynamics. Instead of inequalities, Jarzynski [ll| gave the exact result connecting nonequilibrium 
^T) • quantities with the equilibrium free energies. In the paper [l3j], Clausius-Duhem processes were analyzed in terms 
of averaging over the ensemble of microscopic realizations. Incidentally, the Jarzynski equality and some related 
inequalities are obtained in [l3l ] . As Crooks showed the Jarzynski equality follows from the assumption that the 
system dynamics is Markovian and microscopically reversible. Note that these two conditions are both adopted in 
many computer simulations. In Q, Crooks also derived a related fluctuation theorem of own significance. Quantum 
counterparts of both the results were given by Tasaki [3911 . In the literature, various derivations and reformulations 



of Jarzynski's equality have been presented (see, e.g., |l4 l32l. |3~I1. l38j and references therein). 

Methods of deriving Jarzynski's equality are often based on a dynamical description within an infinitesimal time 
scale. For instance, we can use some kind of master equation [l^, Hij , certain forms of deterministic dynamics 0, H|| , 
or the time- reversal symmetry within a description by time-dependent Hamiltonians @, 0]- However, a matter of the 
Jarzynski equality is in dealing with quantities only at the initial and final points, without explicit reference to the 
passage of time. Hence, we may be interested in obtaining Jarzynski's equality with only discrete state changes. An 
analysis of Jarzynski's equality within a unitary time evolution was given in the papers [2lL [39j . On the other hand, 
reversible unitary transformations form a very special class of possible state changes. The formalism of quantum 
operations is one of basic tools in studying dynamics of open quantum systems. This formalism is especially well 
adapted to describe discrete state changes. Due to such a property, quantum operations are widely used in quantum 
information theory |23j . Deterministic processes are represented by trace-preserving operations also known as quantum 
channels. It is interesting to examine the Jarzynski equality from the viewpoint of quantum-channel approach. 

The aim of this paper is to derive Jarzynski's equality and some related results under assumption that a change of 
system state is described by unital quantum channels. The paper is organized as follows. In Section|TTl the preliminary 
material is reviewed. Basic notions of the quantum operation techniques are recalled. Joint probability distribution for 
measurements statistics is discussed as well. Hence, the averaging rule for considered topics is obtained. In Section ITjTl 
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we get Jarzynski's equality under condition that the considered process is represented by a unital quantum channel. 
The derivation is based on a general formal result, which might be useful in other contexts. Further, we formulate a 
theorem of Tasaki-Crooks type for bistochastic map and its adjoint. Using the developed formalism, we also discuss 
a heat transfer between two quantum systems. Initially, the combined system is assumed to be in the product state of 
two particular densities of a special form. In Section IIV1 some equalities with unital channels in the case of feedback 
control are obtained. First, we examine the case of error-free feedback control. Second, we extend the formulation to 
a feedback control with classical errors. In Section [Vj we conclude the paper with a summary of results. 

II. DEFINITIONS AND NOTATION 

In this section, the required material is presented. First, we briefly recall basic notions of the formalism of quantum 
operations. Second, we describe a general form of averaging procedure with a joint probability distribution. This 
procedure will be assumed in nonequilibrium relations such as Jarzynski's equality and the Tasaki-Crooks fluctuation 
theorem. 

A. States, operators, and quantum channels 

Let C(H) denote the space of linear operators on finite-dimensional Hilbert space T-L. By C s . a .{T-L) and £+(H), we 
respectively mean the real space of Hcrmitian operators and the set of positive ones. For any X G we put 

| X | G £ + (TL) as a unique positive square root of X. The eigenvalues of |X| counted with multiplicities are the 
singular values of X. The spectral norm HX^ is defined as the largest singular value of X. For arbitrary X, Y G C("H), 
we define their Hilbert-Schmidt inner product by [43] 

(X,Y) hs :=Tr(X+Y) . (1) 

We now consider a linear map $ : C(Ha) — > £{Hb) that takes elements of C(%a) to elements of C(%b)- To 
describe a physical process, this map must be completely positive pi. l23j|. Let id_R be the identity map on jC^r), 
where the space Tin is assigned to a reference system. The complete positivity implies that Cg> id^ transforms each 
positive operator into a positive operator again for each dimension of the extended space. Any completely positive 
map can be written in the operator-sum representation. For all X G C(1~La), we have 

$ ( x ) = E M K » XK i ■ ( 2 ) 

Here, the Kraus operators K M map the input space T-La to the output space Kb. To each linear map $ : C(T-La) ~ > 
C(Hb), w e assign its adjoint H : C(Hb) —> £(Ha) by the formula [3[ 

(*(X),Y> h8 = <X,$t(Y)> hs) (3) 

which holds for all X G C(Ha) and Y G £(Hb)- For the linear map its adjoint is represented as 

$t(Y) = £ KtYfV (4) 

In general, the state of an open quantum system is described by a density matrix, which is a positive operator of unit 
trace. The input density matrix p A is mapped to the output <&(p A ) G C+^Hb)- To be consistent with probabilistic 
interpretation, the map $ should obey the condition [23| 

E K PV<^> (5) 

where TLa denotes the identity operator on "Ha- Quantum operations are maps of the form @ under the restriction 
([5|). Deterministic processes are described by trace-preserving operations, for which the inequality ([5]) is saturated 
and, herewith, Ty(&(p a )^) = 1. These maps are usually referred to as quantum channels j2;]. Non-trace-preserving 
quantum operations occur in probabilistic protocols such as unambiguous discrimination [27] or probabilistic cloning 
8]. To obtain the output density matrix for probabilistic operations, we merely rescale the output as 



p B := Tr($(p A )) 'HP a) 



(6) 
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Except for trace-preserving maps, the denominator in ([6]) generally depends on the input pa- In this way, therefore, 
we do not always obtain a trace-preserving operation. On the other hand, the macroscopic dynamics is deterministic. 
For these reasons, we further focus an attention on quantum channels. Let us consider the operator sum 

<p(t A )=j2^i- (7) 

Due to the positivity and the definition of the spectral norm, we can write 

< ||$(1a)||oo 1b ■ (8) 
Assume that the operator ([7]) is a multiple of 1b- For trace-preserving maps, this condition gives 

$(1 A ) = ^ U ■ (9) 

VLB 

Here, the integers (1a = dim^^) and ds — dim(Hs) are dimensionalities of Ha and Hb, respectively. When the 
input and output spaces are of the same dimensionality, the formula ^ gives = 1_b- The latter is commonly 

said that the map is unital Q. For the simplest system, called quantum bit, the depolarizing and phase damping 
channels are both unital, whereas the amplitude damping channel is not [231. N ote that the depolarizing channel, 
which represents a decohering qubit, has interesting entropic characteristics |29j . Unital trace-preserving maps are 
often referred to as bistochastic 2}. For bistochastic quantum channels, the formula © becomes equality with 
||$(1a)||oo = 1. It is easy to check that the adjoint of bistochastic map is bistochastic as well. In the following, we 
derive a useful auxiliary statement about those trace-preserving maps that satisfy the condition (J9)). 



B. Joint probability distribution and averaging rule 



The basic idea of statistical physics is to represent a macroscopic situation of interest by an ensemble of its mi- 
croscopic realizations [30| . The representative ensemble is constructed in such a way that all the realizations arc 
consistent with the actual macroscopic conditions. For a quantum version of Jarzynski's equality, the corresponding 
framework was explicitly developed by Tasaki [39j. Developing course of the paper [24| . we will initially pose an 
approach for arbitrary observables A £ C s . a .^HA) and B e £ s . a .("Hs). Up to substitutions, our formulation follows 
[39(. In terms of the eigenvalues and the corresponding eigenstates, spectral decompositions of the two observables 
are expressed as 

A = V ai \ ai ){ ai \ , (10) 

z *% 

b = £aI^XM • (ii) 

* — '3 

In both the decompositions, the eigenvalues are assumed to be taken according to their multiplicity. In this regard, 
we treat aj and bj as the labels for vectors of the orthonormal bases and Suppose the system evolution 

is represented by quantum channel If the input state is described by eigenstate |oj), then the channel output is 
<&(|aj)(aj|) . Then the probability of being the state \bj) is calculated as 

p(6^H^(h>Ml&i> ■ (12) 

This quantity is the conditional probability of the outcome bj given that the input state was ) . Due to the 
preservation of the trace, we then obtain 

^,p(b j \a i ) = Tr(<S>(\a i )(a i \)) =1. (13) 

Thus, the standard requirement on conditional probabilities is satisfied with any quantum channel. Further, we 
suppose that the input density matrix p A has the form 

Pa = y^.p{ai)Wt){a l \ , (14) 

where J2iP( a i) = 1- Of course, the operator (fl"4|) could be rewritten as a function of the observable A. According to 
the Bayes rule, one defines the joint probability distribution with elements 



p(ai,bj) = p{ai) p{bj\ai) 



(15) 
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This is the probability that we find the system in z-th eigenstate of A at the input and in j-th eigenstate of B at the 
output. Consider a function /(a, 6) of two eigenvalues. Extending Tasaki's approach [39| . we define an average 

(f(a,b)) : X P(<*i A) A) ■ (16) 

^3 

Here, angular brackets in the left-hand side signify averaging over the ensemble of possible pairs of measurement 
outcomes. Between these brackets, we will usually omit labels of the involved variables. In Section IIV1 however, we 
will consider enough complicated protocols. There, the labels will be all indicated for clearness. 

In general, the average (|16[) does not correspond to quantum-mechanical expectation values. To emphasize the 
point, Tasaki j3(| gave the term "classical average" for (jTHJ) . Namely, the right-hand side of Eq. (|16l) corresponds 
to specific physical meaning. In two simplest cases, however, the average (|16p concurs with the quantum-mechanical 
expectation value. Let g(a) be some well-defined function. Due to (TT5)) and (H~4l . we directly obtain 

(5(a)) = J2.9{ai)p(ai) = ^^^PA^i) = Tr( ff (A)p A ) . (17) 
Using the map linearity and the formulas (fT2|) and (fT5|) . we also write 

(9(b)) =X>fe)£ .J>(«iWil*(l«iX«il)l&i> = n(<7(B)$(p A )) • (18) 

3 1 

These conclusions are essential for a physical interpretation of all the further calculations. They should be kept in 
mind in the following. 



III. RELATIONS OF JARZYNSKI AND TASAKI-CROOKS TYPES 



In this section, we obtain some results connected with Jarzynski's equality and the Tasaki-Crooks fluctuation 
theorem. The derivation will mainly be based on the condition that the representing quantum channel is unital. 
Further, we consider a heat transfer between two quantum systems, when state changes of the combined system are 
also governed by unital quantum channels. 



A. Jarzynski's equality with unital quantum channels 



Before obtaining Jarzynski's equality, we will formulate a mathematical result in more abstract form. Let us consider 
the case, in which the input density matrix is expressed as 

:= Tr(e-<* A ;f V aA . (19) 

Functional form of such a kind is related to the state of thermal equilibrium in the Gibbs canonical ensemble. The 
following exact relation could be applied beyond the context of Jarzynsky's equality. 

Proposition 1 Let A £ <Cs.<i.(%a)> B G Cs.aX'H.B), md let a and f3 be real numbers. Suppose the input state is 
described by density matrix h!9\) . If the quantum channel $ : C(iHa) C^Hb) satisfies the condition {P|) ; then 

<«,(<—*)) = <»> 

Proof. Using the linearity of the map $ and the condition (|9"]l. we obtain 

j^p&ioo = (b 3 \^2 % H\^)(^\)\b 3 ) = (bjiHum = ^ . (21) 

Taking p(a,) = Tr(e~ aA ) 1 exp(— acij) in (TT5|) . we then express the left-hand side of (|2U)) as 

^l^Sf P(>il«0 «p(°*-/%) = ^ry^y I>pH%) f B ■ (22) 
The latter term concurs with the right-hand side of (|2U| . ■ 



5 



The condition (j9)) has allowed to evaluate the sum (|2Tj) in a closed form. Hence, the claim ((20)) immediately follows. 
Basing on the result (I2TJ1) . the Jarzynski equality can easily be derived. Concerning the sum (|2ip . a comment is 
required. This formula should not be confused with the constraint (fT3| . which takes the sum with respect to j at 
fixed i. The distinction may be illustrated with the equiprobable distribution p(a,i,bj) = (fiA^s) -1 ■ Hence, we have 
the conditional probabilities p(bj\a,i) = dg 1 and, herewith, the right-hand side of (|21j) . In the context of quantum 
channels, this example can be realized as follows. Assuming <1a < ds, we take an isometry V : Ha — > Hb such that 
the orthonormal set |V|aj}} is mutually unbiased with the basis In other words, one gives |(&j|V|ai)| = d^ 2 . 

The trace-preserving map $ is then defined by ([2]) with a unique Kraus operator K = V. In many cases of interest, 
the system dimensionality is not altered during a physical process, i.e. d& — d B . 

Let us proceed to a physical situation corresponding to the context of Jarzynsky's equality [l5l |37| . We assume 
that a thermally contacted system is acted upon by an external agent. This agent operates according to a specified 
protocol. Hence, the Hamiltonian of the system is time-dependent. The principal system is initially prepared in the 
state of thermal equilibrium with a heat reservoir. Following Tasaki [39j, we will firstly assume that the reservoir 
temperature is also dependent on the time. The parameters /?o and (3\ give the inverse temperature of the reservoir 
at the initial and final moments, respectively. Thus, the initial density matrix is 

w„(A>) = Z (/3 )- 1 e-' 3oHo , (23) 

in terms of the initial Hamiltonian Ho and the corresponding partition function Zq(Po) — Tr(e _,3oHo ). We further 
suppose that the transformation of states of the system is represented by quantum channel $ with the same input 
and output Hilbert space. In general, the final density matrix $(u>o(/3o)) wm enough differ from the matrix 

Wl (/3 1 ) = Z 1 (/3 1 )" 1 e- /3lHl , (24) 

corresponding to equilibrium at the final moment. Here, the partition function Zi(fti) — Tr(e~^ lHl ) is expressed in 

terms of the final Hamiltonian Hi. By {ffm } ancl { £ n }i we respectively denote eigenvalues of the Hamiltonians Ho 
and Hi. Taking dA = ds implies that the map is unital. By obvious substitutions, the formula (|20j) then gives 

(exp( fte <«>- A£ <»)> = |||, (25) 

provided that the channel $ is unital. Assuming unitary evolution, the relation (|25[) has been derived by Tasaki 39] . 
So, we have extended an important formulation to unital quantum channels. 

In the considered context, the term w nm = e n — £m is treated as an external work performed on the principal 
system during a process. At constant temperature, i.e. when /3q = Pi — (3, we therefore have 

^exp(-/3w)^> = cxp(-/3AF) , (26) 

since Ft(/3) — — /3 _1 \nZt(f3) for t = 0,1. The result ([26)) relates, on average, a non-equilibrium external work with 
the difference AF = F\ — Fq between the equilibrium free energies. This formula is the original Jarzynski equality 

[HIS- As 

a consequence, the basic inequality of thermodynamics can be obtained. Combining (f26| with Jensen's 
inequality for convex function x i— > exp(— flx) leads to 

exp(-0(w)) < exp(-fiAF) . (27) 

As the function x H> exp(— fix) decreases with x for positive (3, the formula (|2T|) gives (w) > AF. Namely, total 
external work will, on average, exceed the difference between values of the equilibrium free energy at the final and 
initial moments. Basing on the Jarzynski equality, Tasaki also discussed some inequalities for the von Neumann 
entropy [39| . The Jarzynski equality has been generalized to the case, when nonequilibrium classical systems are the 
subject of a feedback control [32|, Hj|. In quantum setting, corresponding version was given in [2~j| . 



B. Theorem of Tasaki— Crooks type for bistochastic map and its adjoint 

Fluctuation theorems are traditionally used in studying stochastic processes. Recent progress in statistical me- 
chanics has lead to a renewed interest in this topic [HI O, [13, H3] • Fluctuation theorems can be used for deriving 
information-theoretic results such as Holevo's bound [17| . For the Tasaki-Crooks fluctuation theorem, a development 
with unital quantum channels is motivated as follows. For trace-preserving map $, its adjoint $t [ s un ital. To make 
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trace-preserving, the map $ itself should be unital as well. For this reason, we focus an attention on bistochastic 
maps, i.e. on unital quantum channels. By %, we denote the Hilbert space assigned to the principal system. Similarly 
to (fT9|) . we introduce the density matrix 

£ B (a) Tr(e- QB rV QB . (28) 

Consider two processes obtained by applying the channel $ to the input g A and the channel $t to the input g B . In 
each of the processes, we can ask for a probability that the difference (aj — bj) takes a certain value. It turns out that 
the two corresponding probabilities obey some relation. For a difference between eigenvalues of the Hamiltonians H 
and Hi, results of such a kind are usually referred to as the Tasaki-Crooks theorem 0,H3- It will be convenient, 
however, to pose a statement in more abstract form. Let us put the notation. Imposing the restriction bj — at = A, 
we fix a difference between eigenvalues of the observables B and A. By P(b — a = A| $, g A ) , we denote the probability 
of this event given that the channel $ represents an evolution of the input g A . We have the following result. 

Proposition 2 Let the quantum channel $ : £(H) — > £-{%) be unital. For all real A, the corresponding probabilities 
satisfy 

e- aA Tr(e- aA )P(b-a = A\<5>,Q A )=Tr(e- aB )P(a-b=-A\<I>\Q B ) , (29) 

where the density matrices g A and g B are respectively defined by H9\) and i28\ ) . 

Proof. Assuming action of the channel $t, we put the conditional probability of outcome given that the input 
state was \bj). Similarly to (|12l) . this probability is written as 

q(a i \bj) = (a i \&(\b j )(bj\)\a i ) . (30) 

The derivation of (|29|) is based on a simple observation. Due to the representations ([2]) and ((4]), we have 

(b|K»(a|Kt|6) = (a|Kt|&)(&|K» , {6|$(|o)<o|)|6) = (a|*t(| 6 )(6|)| a > , (31) 

with arbitrary \a), \b) S H. For all i and j, therefore, one satisfies 

p(bj\ai) = q(ai\bj) . (32) 

The matrix (TT91 has eigenvalues p(a.i) — Tr(e~ aA ) 1 exp(— aai). Let A/"(A) be the set of ordered pairs (j, i) such that 
bj — cii = A. By (|32|) . the left-hand side of (|29|) is then rewritten as 

} J exp(-a6 i + aai)Tr(e~ aA ) p(ai)p(b 3 \ai) = exp(-abj)p(bj\ai) 

(j,i)eAf(A) (j,i) eA/" (A) 

= Tr(e-« B ) l(bj)q( ai \bj) . (33) 

(j.i)eAT(A) 

Here the numbers q(bj) = Tr(e _aB ) 1 exp(— abj) are eigenvalues of g B . For all the ordered pairs (j, i) 6 A/"(A), we 
have ai — bj = — A. Combining this with the definition (l3Tfl) shows that the right-hand side of (|3"3")) concurs with the 
right-hand side of (US]). ■ 

In line with the above derivation, we can also write some relation with the completely mixed state p f = l/<i, where 
d = dim(H). When p(ai) — q(bj) = 1/d, the formula (|33)l is replaced with 

p{a i )p{b J \a l ) = Y l( b ]) q(a t \bj) ■ (34) 
OVt)eAf(A) <j,i)eM(A) 

Hence, we have the corresponding relation 

P(b-a = A|$, p,) =P( a -b = -A|$+, p,) . (35) 

It seems that the matrices of the form (|19[) and the completely mixed matrix p^ are the only two forms, for which a 
closed relation between the two probabilities can be written. 

The statement of Proposition^ is expressed in a general form. It immediately leads to a theorem of Tasaki-Crooks 
type. Let us consider a thermally insulated system acted upon by an external field. As above, the operators Ho and 
Hi are respectively the initial and final Hamiltonians. For a unital quantum channel $, we apply this channel itself to 
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the equilibrium state ((23)) and the adjoint $^ to the equilibrium state (|24|) . Differences w nm = el 1 ' — are treated 
as possible values of an external work performed on the system during the former process. By obvious substitution, 
the formula (|2T)|) gives 

e-P w Z (j3)P(sW =w\*,u Q ) =Z 1 ()3)P(s^ -e^ = -u^uii) , (36) 

with the inverse temperature /3 of heat reservoir. Using AF = Fi — Fq, the relation (|36j) can be rewritten in the form 



P( £ (i)_ £ (0) =ti ,|$ jfa , ) 
p( e (0) _ £ (1) ~ _ w |$t ;U ,^ 



= exp(/3w - j3AF) , (37) 



when probabilities are non-zero for taken w. If the channel 4> represents a unitary evolution, then its adjoint $^ 
represents the inverse unitary evolution. In the case of unitary transformations, the exact relation (|36[) with the two 
probabilities has been derived by Tasaki [3{|. It is a quantum analog of previous Crooks' formulation @. In the 
literature, the above statement is often referred to as the Tasaki-Crooks fluctuation theorem 0, [Hj]. Thus, we have 
obtained an extension of the Tasaki-Crooks fluctuation theorem to unital quantum channels. Of course, the Jarzynski 
equality is closely related to the Tasaki-Crooks theorem. Since this issue is well discussed in the literature, we refrain 
from presenting the details here. 



C. Notes on heat transfer between two systems 

We now apply the above results to a heat transfer between two quantum systems. The following analysis is an 
extension of related results of the paper (39j . It is convenient, however, to put a derivation in a more abstract manner. 
Consider two systems with the Hilbert spaces T-La and %b, respectively. So, we write the initial density matrix of the 
combined system as 

Qab ~ Qa{v) ® Q B (0) = TV(e-« A )- 1 Tr(e-^ B )- 1 ® e"^ B , (38) 

with A £ C s .a.(HA) and B £ C s . a .(T~tB)- Here, we used the matrices (fT§|) and (|2"5)l . but the latter with /? instead of a. 
We also assume that state changes of the combined system are represented by the channel with the input space 
T~i-AB = Ha ®T~Lb and the output space He- Let us take quantum channels that satisfy condition of the form (JSJ, 
namely 

a,(t AB ) = ^ii C) (39) 

where \ab = 1a ® 1b and lc are the corresponding identities. The operator aA Cg> 1^ + 1a ® /3B is Hcrmitian for 
real a and /3. It has eigenvalues ua,i + (3bj and eigenstates \a,ibj) — \ai) ® \bj). In terms of this operator, the matrix 
(J35J reads 

q ab :=Tr^exp(-o;A(8)lfl-lA(8) / 0B)^ exp(-aA ® 1 B - 1a ® /SB) , (40) 

since the summands aA (g> 1b and 1a <8> /SB commute. For arbitrary C € £ s .a.(%c), the following conclusion can be 
written as a variety of (|20|) . If the quantum channel ^ : C(H A b) — > £(Hc) satisfy ([39]) and the input is given by 
I, then 

exp(aa + &b - c)\ = ? , . (41) 

V ' I d c Tr(e-« A ) Tr(e^ B ) V ' 

In the case He — Hab and C = a A <8>1b + 1a<8>/3B, the formula (j4]~|) becomes 

^exp(a(a - a') + 0(6 -&'))) = 1 . ( 42 ) 

where the set {aa- + (3b'A denotes the spectrum of C. The condition ([39)) is reduced here to ^(Iab) — Iab- 

We now take two systems initially prepared in the equilibrium at the inverse temperatures /3q and ft , respectively. 
Their density matrices are therefore written as (l23l) and (|24"]) . whence the product u>o(ft) ® u>i(ft) is the input total 
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state. Suppose that the systems further interact via unital quantum channel 'J. By obvious substitutions into (|42l) . 
for the described process we obtain 

(exp(/3 (e (0) -e' (0) )+/3i(£ (1) = 1 . (43) 

Assuming unitary evolution, this result has been given in j39j. Thus, we have extended the previous result to unital 
quantum channels. In the paper [39| . the relation (|43[) is presented with a certain physical interpretation. Let us 
consider the quantity 

AS := (/3 (e' (0) - £ (0) ) + - . (44) 

The terms (e'W — e^ * 1 ) and (e'^ 1 - 1 — eW) give, on average, the change of self-energy of the corresponding systems. 
Suppose that changes in the inverse temperatures of the systems are sufficiently small and a contribution of interaction 
energy to the total entropy is negligible. In such a situation, the quantiy AS* estimates an averaged change of the 
total entropy [39]. Combining (1431) with Jensen's inequality for convex function x i-4 exp(— x), we obtain 

exp(-AS) < 1 , AS > . (45) 

Thus, we have arrived at the well-known inequality of thermodynamics. In the paper |39| . Tasaki also discussed a 
way for constructing more detailed bounds on AS* from below. We only emphasize here that the consideration can 
quite be proceeded with unital quantum channels. 

IV. JARZYNSKI-SAGAWA-UEDA RELATIONS WITH UNITAL QUANTUM CHANNELS 

In this section, we develop some of the above results in the case, when the agent makes measurement followed by a 
feedback. First, error-free feedback control is considered. Second, we analyze the case, in which classical errors occur 
in the measurement process. 

A. Error-free feedback control 

We assume that the agent performs a quantum measurement and further changes the protocol according to the 
measurement outcome. For classical systems, this topic has been considered b y S agawa and Ueda [32, [33j. For 
quantum systems, relations of such a kind were examined by Morikuni and Tasaki [21] . Let us recall required material 
on quantum measurements. In general, the quantum measurement is described by a set {N^} of measurement 
operators, acting on the space of the measured system [23]. If the pre- measurement state is described by p, then the 
probability of /x-th outcome is Tr(N^N^p). The corresponding post-measurement state is described by density matrix 

TxfNtN^-XpNt . (46) 

It must be stressed that number of measurement outcomes can exceed dimensionality of the Hilbert space. For instance, 
such a possibility is crucial for unambiguous discrimination of nonorthogonal states (27j . The set of measurement 
operators should satisfy the completeness relation 

E N X = 1- (47) 

When the measurement operators are mutually orthogonal projectors, the above scheme obviously leads to traditional 
projective measurements. We will also say about quantum measurements that satisfy the condition 

52^1 = 1. (48) 

If the measurement operators obey (|48p. then the corresponding measurement will be said to be unital. This term 
agrees with the definition of unital maps, for which the right-hand side of ([7]) becomes the identity operator. Of 
course, all the projective measurements satisfy the condition (|4"8")l . Moreover, for Hermitian measurement operators 
this condition concurs with the completeness relation (|47[) . 

Let us pose formally the protocol with free-error feedback control. For the sake of simplicity, we will assume the 
same input and output space for all adopted operations. Further, we put an appropriate number of observables 
B M 6 C s .a.{'H) 1 each with its eigenbasis )}. We also define the corresponding density matrices Q^{a) by the 

formula (|28p with instead of B. We shall consider the following procedure. 
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(i) At the first stage, the agent applies quantum channel $ to the input g A given by (IT9l) . 

(ii) At the second stage, one performs a quantum measurement on the output Q(g A ). For /i-th outcome, its 
probability is p{p) = Tr(l\|t N A1 $(g A )) and the post-measurement state is p([i)~ 1 N^$(g A )Nj i . 

(iii) At the third stage, the agent applies the prescribed quantum channel 4" „ to /i-th post-measurement state given 
that /i-th outcome occurs. 

(iv) At the fourth stage, one measures the observable B M on the third-stage output p(/i) _1 5' Al (N A1 , I>(g yl )Nj 1 ). With 
this pre-measurement state, the outcome bj^ is obtained with the probability 

p( M )- 1 (^ ) |^(N Al $( eA )Nt)|^ ) ) . (49) 

Multiplying (|49|) by p(fi), i.e. by /i-th outcome probability, we obtain the probability of the outcome bj for the 
input (|19l) . The latter probability can be represented as the sum 

]T pK,^) = J2Mp( b f\ a i) ■ (so) 

Here, we introduce the conditional probability 

p(bf\ ai ) = (6^ ) |^(N A1 $(|a 4 )(a 4 |)Nt)|^) , (51) 
given that the input was \ai). These probabilities satisfy the required condition 

^^ ) M=^Tr(NtN A1 $(|a 4 )(a J |))=l . (52) 

Here, we have used the preservation of the trace by channels 4^ and $, the cyclic property of the trace and the 
completeness relation (|47l) . For the considered scenario, our main result is formulated as follows. For clearness, all 
the labels of involved variables will be explicitly indicated between the angular brackets. 

Proposition 3 Let the above protocol be applied to the input M9\) , and let the quantum channel $ : C{%) — > £■{%) 
be unital. For arbitrary quantum measurement {N^} and quantum channels 4*^, we have 

(■^Sr) e X p(^-^ ) ))=E M ^(^ ) («)^( N . N t)) ■ (53) 
Proof. Since the channel $ is unital and (b^\ exp(— ctb^) = (b^\ cxp(— aB^), we first write 

^piaMbfh) expK-a^>) = ^[g^ (bf ) \bf) 

= (6f |^(a)^(N M Nt)|6f ) . (54) 

Summing the right-hand side of (f54| with respect to j, we obtain Tr(g^ (a) ^(N^Nj,)^. The latter leads to (j53|) 

after further summing with respect to \x. I 

The statement of Proposition [3] is written in a general form. We now apply this result to a thermally insulated 
system, which is operated by the agent with feedback control. For each of the possible ways, we introduce the 
corresponding Hamiltonian H M and equilibrium state 

uM = Z^)- l e-^ , ZM = Tr(e-^) . (55) 

The protocol is applied to the input u>o(/3), whereas the averaging is taken over the final states u}^{f3) with /i ^ 0. 
Let us introduce the associated work Wftm = — Em ■ By obvious substitution, we rewrite the equality (|53l) as 

(exp(-/3^ + - F ))) = 7 , (56) 
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where F^(f3) = —(3 1 lnZ^/?) and the parameter 7 is defined by 

7:=£^(«mG9)*/.(N m N],)) • (57) 

The formula (|56[) is a quantum counterpart of one of the results originally given by Sagawa and Ueda [32j . Assuming 
unitary transformations of quantum states, the relation (|56|) was derived in [2l| . For two-time projective measurements 
and any trace-preserving map between them, an analog of (I56|) was formulated and applied in [l7j . Thus, we have 
generalized an important nonequilibrium equality to the case, when the unital channel $ acts at the stage (i) and 
arbitrary channels ^ act at the stage (iii) . Take the case, when the same channel *f> is applied for all the measurement 
outcomes and final states are always compared with u>i(/3). That is, the agent makes a measurement, but does not 
perform any feedback. If the channel ^ and the measurement are both unital, then we have 7 = Tr(a>i(/3)) = 1. 
Namely, we have arrived at the Jarzynski equality (|26[) with an intermediate quantum measurement. Some physical 
consequences of ([56| are discussed in [2l[. This result allows an extension to feedback control with classical errors. 



B. Feedback control with classical error 



The above result can be extended to the case, when measurement outcomes are registered with a certain amount 
of randomness. Assuming purely classical nature of errors at the stage (ii), it is represented by probabilities r(y\[x). 
The quantity r{y\ji) is the conditional probability of mis-interpretation of actual /i-th outcome as registered ^-th one. 
For the following formulations, some additional notions are required. Recall the concept of mutual information. The 
pointwise mutual information is defined as [9j 

I»u ■= In . (58 

It can take positive as well as negative values, vanishing for z/) = p(fi)p(v). Averaging (|58l) with the joint 
probability distribution, we obtain the mutual information 

<V>=£K/^)m^. (59) 

Due to the Bayes rule, the pointwise mutual information can be rewritten as 

r(jt\v) r{v\i±) 

I^ v =m — — - =ln — —— . (60) 
PW P{v) 

This measure is obviously symmetric. Other properties of the mutual information and its relations to entropic measures 
are extensively treated in information theory books such as @. 

In the case of feedback control with classical errors, we should average with the joint probability distribution 

p(ai,H,b^) = p(a i )p(fj,,b < f ) \ai) . ( 61 ) 
in which the latter conditional probability is written as 

pGMf |a0 = r(HM)(^ ) |^(N^(|a i )(a i |)Nt)|6$ !/) ) . (62) 

For error-free feedback control, we have r(u\fi) — <5„ M . Here, the right-hand side of (|62p is non-zero only for v = ji, 
when it is reduced to the right-hand side of (fBTj) . Similarly to (|52"|) . we immediately obtain 

p(»,bf ) \a i ) = l, (63) 

in view of r{y\fi) = 1. We can now formulate two nonequilibrium equalities in the presence of feedback control. 

Proposition 4 Let the above protocol be applied to the input H9\) , and let classical errors at the stage (ii) be repre- 
sented by the conditional probability r{y\p,). If the quantum channel $ : C(l-L) —> C(T-L) is unital, then 



Tile 



-aA\ 



Trie 



-qB„ 



exp( aai -abf)) =^r(z,| A1 )Tr(^ ) (a)^(N M Nt)) . (64) 
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If the quantum channels : C{T-L) — > C(H) and the quantum measurement {N M } are also all unital, then 

( Tr ( e _ aB „) exp(aa, - abf - I Ufl ) J = 1 . (65) 
Proof. Similarly to the proof of Proposition |3l we first obtain the relation 

^Je?) E^M^&fk) exp(a ai -a6^) =r(H/i)Tr(^ ) (a)^(N,Nt)) , (66) 

provided that the channel $ is unital. After summing with respect to and v, we get the first claim (1641) . Multiplying 
(I66p by exp(— 7,,^) = p(z/)/r(^|^i) and further summing with respect to /i, one yields 



^I-ob!) Ey^fC *)^' 6 ! 1°*) expK -a&f -7^) = 



(67) 



under the condition (|48|) . If = 1 for all is, summing of (|67|) with respect to v finally gives J2i/P( u ) = L ■ 

It is essential for and (f6"5"j) that corresponding quantum operations be unital. The proof of the result 
assumes this property only for the channel at the first stage (i). The proof of the result (|6"5")) has also assumed this 
property for the measurement at the stage (ii) and for all the channels at the stage (iii) . Let us proceed to a thermally 
insulated system, which is operated by the agent with feedback control. In the notation of this section, the formula 
gives 

(exp(-/3w^ + P{F V -F )))=j, (68) 

where the parameter 7 is written as 

^ : =E^ r ^i^ Tr (^ (/3)VI/ ^ N ^ N ^) • (69) 

Replacing 7 with 7, the result (|68|) is completely similar to (|56l) . For error- free feedback control, the term 7 concurs 
with 7 due to r(v\n) — 5^^. The formulas (1551) and (|6"5|) express the parameters in a closed form. Note that these 
parameters are experimentally measurable quantities. The original Sagawa-Ueda formulation of Jarzynski's equality 
with feedback control has been tested experimentally [4(| ■ The results (j51)f and (|68p have been proved under the same 
condition that the channel $ is unital. When the quantum operations at the stages (ii) and (iii) are unital as well, 
we further have 

Pw^l + /3(F U -F )-I vll ))=l. (70) 



exp( 

This is a quantum counterpart of one of the results obtained by Sagawa and Ueda [32| . In the paper (2l| , the 
relation (|70[) has been presented within a unitary evolution under assumption that the measurement obeys 



The formula (|68[) has also been given in 2l| with unitary transformations of the states. Thus, we have extended 
quantum Jarzynski-Sagawa-Ueda relations to unital quantum channels. Using the equalities (|56p . (I68[) . and (|70l) . 
some interesting inequalities and observations can be obtained. We refrain from presenting the details and refer to 



section 4 of 21(. One of important issues is the maximal work that can be extracted from multi-heat baths with 
the assistance of discrete quantum feedback control [3l|. In the paper [3l|, the maximal work is characterized by a 
generalized mutual information, which is different from the above mutual information. 



V. CONCLUSIONS 



Recent advances in studying nonequilibrium processes are connected with exact relations such as the Jarzynski 
equality and fluctuation theorems. We have considered this issue from the viewpoint of quantum operation techniques. 
Any details of the passage of time are not requested. In the paper, we have emphasized the following point. It is 
essential that representing quantum channels be unital. Assuming this property, the Jarzynski equality is still valid. 
A fluctuation theorem of Tasaki-Crooks type has been formulated for bistochastic map and its adjoint. With unital 
quantum channels, we also apply the formalism to the problem of heat transfer between two quantum systems. Some 
equalities with unital channels have also been derived in the case of feedback control. Error-free feedback and feedback 
with classical errors are both considered. Hence, quantum Jarzynski-Sagawa-Ueda relations have been generalized 
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to unital quantum channels. The obtained expressions may be useful in experimental tests. In general, the presented 
consideration can be viewed as a supplement and development of the previously given facts about nonequilibrium 
equalities. It could provide a novel insight into a dynamical behavior required for the validity of Jarzynski's equality 
and related results. 

Note added 

After the present paper was submitted I learned about the recent work [l(, in which the significance of unitality 
is emphasized as well. In the work [lj], the authors formulate a general fluctuation theorem and further show that 
some previous results follow from this theorem. Our formulations and derivation methods are different and, in certain 
respects, complementary to that is given in pj. The authors of [l[ also describe results of a related experiment with 
superconducting flux qubits. 
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